Piezoelectrics by Design: A Route through Short-period Perovskite Superlattices 
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Using first-principles density functional theory, we study piezoelectricity in short-period super- 
lattices made with combination of ferroelectric and paraelectric components and exhibiting polar 
discontinuities. We show that piezoelectric response of such a superlattice can be tuned both in 
terms of sign and magnitude with a choice of its components. As these superlattices with non- 
switchable polarization do not undergo ferroelectric transitions, we predict them to exhibit a robust 
piezoelectric response with weaker temperature dependence compared to their bulk counterparts. 
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With advances in experimental techniques, it is now 
possible to create superlattices of various of oxide per- 
ovskites (general formula ABO3) through layer by layer 
epitaxial growth. Such superlattices have attracted lot of 
attention in recent years due to the presence of unusual 
properties .'^^"'^ In particular, these superlattices may ex- 
hibit polar discontinuities leading to the so-called "polar 
catastrophe" , a topic of much discussion in recent time 
in the context of LaAlOa/SrTiOa thin films.^^^V^ 

Superlattices with polar discontinuity were studied by 
Murray and VanderbillP from first-principles for super- 
lattices consisting of ferroelectric (FE) and paraelectric 
(PE) components, and the "polar catastrophe" was found 
to to result in response of system in terms of polar dis- 
tortions at the interface. This was extended in our recent 
workpSlto superlattices with ultra-thin (1 unit-cell) com- 
ponents of FE and PE layers. Our studyl^l carried out 
on a series of such superlattices constructed from differ- 
ent PE and FE components showed that even for ultra- 
thin superlattices, systems containing polar discontinu- 
ities have strongly broken inversion symmetry and large 
non-switchable polarization. The magnitudes of the po- 
lar distortions are as large as those in the corresponding 
bulk FE materials. Evidently, this class of superlattices 
form promising new material with interesting properties. 

In a recent study, tunability of piezoelectric proper- 
ties of short period ferroelectric superlattices with strain 
was investigated.'HlHere we demonstrate tunability of ro- 
bust piezoelectric properties for superlattices formed of 
FE and PE layers through (i) polar discontinuity and (ii) 
chemical control, using first-principles density functional 
theory (DFT) based calculations. 

We consider two representative systems: 
LaAlOa/PbTiOa (LAO/PTO) and KSbOa/PbTiOg 
(KSO/PTO) and design composite super- 
lattices formed out of these two superlat- 
tices, such as KSbOa/PbTiOa/LaAlOg/PbTiOa 
(KSO/PTO/LAO/PTO). The first two systems are 
composed of ahernating III-III (A+'''B+303)/II-IV 
(A+^B+'^Oa) and I-V (A+iB+^03)/II-IV perovskite 
layers respectively, while the last system is composed of 
alternating TV/II-IV/III-III/II-IV sequence. 

Our calculations are carried out in the plane wave 



basis within the framework of local density approxima- 
tion (LDA) of DFT using the Vienna Ab-initio simu- 
lation package (VASP).I2^H^ The positions of the ions 
are relaxed towards equilibrium until the Hellmann- 
Feynman forces become less than 0.001 eV/A. A 6x6x4 
Monkhorst-Pack k-point mesh and 450 eV plane-wave 
cut off are used for the calculations. The polarization 
is calculated using the Berry phase approach^^ of mod- 
ern theory of polarization. To determine piezoelectric re- 
sponse, strains of ± 1% are applied. The superlattices 
are constructed along the crystallographic z direction, 
and the piezoelectric coefficient relevant to applications 
is Bzz, strain is therefore applied along the z direction. 

Fig.l shows equilibrium geometries obtained by relax- 
ation of the atomic positions and the out of plane lat- 
tice constants (c), with the in-plane lattice constant (a) 
fixed to the LDA bulk lattice constant of PbTiOs in the 
tetragonal phase. The presence of polar discontinuities 
drive both LAO/PTO and KSO/PTO strongly polar, 
as found in RepSl xhe opposite orientation of the po- 
lar discontinuity in case of LAO/PTO and KSO/PTO 
gives rise to off-centric movements of Pb and Ti that are 
oppositely oriented in two cases. The off-centric move- 
ments driven entirely by the polar discontinuity are non- 
switchable. The third superlattice built up with combi- 
nation of LAO/PTO and KSO/PTO superlattices, im- 
pose oppositely directed polar fields on the PTO layer, 
as shown in Fig.l. This leads to the interesting situa- 
tion, where the off-centric movements of Pb and Ti ions 
become negligibly small and the cations in the PE lay- 
ers move significantly with the movements of La and Al 
being opposite in direction to that of K and Sb. The off- 
centric movements of La, Al, K and Sb measured with 
respect to the centers of their oxygen cages are found to 
be as large as 0.25 A, 0.16 A, -0.21 A and -0.17 A respec- 
tively. This may be compared with the large displace- 
ments of Pb and Ti in case of LAO/PTO and KSO/PTO, 
which are found to be -0.25 A and -0.15 A respectively 
for LAO/PTO and 0.29 A and 0.18 A respectively for 
KSO/PTO. This gives rise to a calculated polarization 
value of 12.79 jjlC/cw? compared to -50.87 ^jlC/cw?' for 
LAO/PTO, -f 57.74 [iCjcrr? for KSO/PTO, reduced by 
a factor of 4 - 4.5 in composite superlattice compared to 






FIG. 1: (Color Online) Equilibrium structures of LAO/PTO, KTO/PTO and the composite superlattice KAO/PTO/LAO/PTO 
(from left to right). The arrows mark the various off-centric movements of the cations measured from the centers of their 
respective oxygen cages. For clarity, only the off-centric movements larger than 0.05 A are shown. The polar discontinuities 
arising in each superlattices are shown in terms of formal charges (p) . 



its components. 

To determine piezoelectric coefficients of the con- 
structed superlattices, in the next step we calculate po- 
larization at unstrained structure (P'''^-'') and at several 
strained structures (P'^), with relaxed atomic positions. 
The slope of the (P*^ - P'^'^f) versus strain curve yields the 
piezoelectric constants, Bzz- Fig. 2 shows such curves for 
LAO/PTO, KSO/PTO and KSO/PTO/LAO/PTO su- 
perlattices. We find the piezoelectric coefficient to be 
large and negative (-2.58) for LAO/PTO, large and pos- 
itive (+4.08) for KSO/PTO and a small value (0.48) 
for KSO/PTO/LAO/PTO. Using combination of vari- 
ous perovskite materials, we are therefore able to design 
piezoelectrics with Czz of opposite sign and varied mag- 
nitudes. 

To have large piezoelectric effect, the system must con- 
tain ions with large Born effective charges (Z*) and they 
should easily move as a result of lattice strain. A piezo- 
electric coefficient can be separated into two parts: a 
clamped-ion contribution evaluated at vanishing values 
of the movements of the ions (u) , and a term that is due 
to the displacements of differently charge sublatticea^^ 
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where fl is the volume, Qe is the electronic charge 
and the subscript k corresponds to the atomic sub- 
lattices. Table -I lists the various contributions for 
LAO/PTO, KSO/PTO and KSO/PTO/LAO/PTO su- 
perlattices. We find large sensitivity of atomic move- 
ments upon strain for Pb, 01 (oxygen at Pb plane) for 
LAO/PTO system, for Pb, Ti, 01 (oxygen at Pb plane), 
05/06 (oxygens at the Sb plane) for KSO/PTO sys- 
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FIG. 2: (P'=-P'''=^) for LAO/PTO (circles) and KSO/PTO 
(squares) plotted as function of applied strain along z direc- 
tion. Inset shows the same plot for KAO/PTO/LAO/PTO. 
The slope of each curve give the measure of piezoelectric co- 
efficient (see text). Note the change in y-scale between the 
inset and the main plot. 



tem. The composite system on the other hand, is far 
more rigid, with du/de values being order of magni- 
tude smaller. Born effective charges, however, change 
only slightly between the composite superlattice and the 



TABLE I: Internal displacement gradients as a function of strain {duh,z/dez), Born effective charges (Z*^) for different ions in 
LAO/PTO and KSO/PTO sup erlatt ices. The numbers in the bracket indicate the corresponding numbers for the composite 
superlattice KSO/PTO/LAO/PTO. Clamped ion contribution to total piezoelectric coefficient (e"), contribution due to internal 
microscopic strain (e^), the total contribution c"+e^ (e™*'^'^'^*) and the piezoelectric coefficient obtained from the slopes in 
Fig.2 {ef^^"^) are also listed. 



LAO/PTO 


KSO/PTO 


k 


duk,z 


Kz 


7* * t)Uk,z 


k 


duk,z 


Kz 


Kz * 'l"^' 

dez 


Pb 


-0.f5 (0.03) 


3.42 (3.05) 


-0.5f (0.08) 


Pb 


0.f4 (0.03) 


2.88 (3.08) 


0.40 (0.09) 


La 


-0.06 (0.00) 


4.7i (3.70) 


-0.28 (0.00) 


K 


O.fO (0.02) 


f.23 (f.35) 


0.f2 (0.03) 


Ti 


-0.08 (0.02) 


5.56 (4.82) 


-0.44 (0.08) 


Ti 


0.f5 (-0.02) 


5.38 (5.52) 


0.80 (-0.f3) 


Al 


-0.04 (-0.06) 


3.69 (3.63) 


-0.f4 (-0.20) 


Sb 


0.f2 (0.05) 


6.04 (5.85) 


0.72 (0.3f) 


Of 


O.fS (0.05) 


-3.73 (-3.39) 


-0.48 (-O.fS) 


Of 


-0.f4 (-0.03) 


-3.ff (-3.43) 


0.44 (O.ff) 


02/03 


0.09 (0.00) 


-2. if (-1.74) 


-0.f9 (0.00) 


02/03 


-0.09 (-0.02) 


-f.6f (-f.68) 


0.f5 (0.04) 


04 


0.04 (0.06) 


-3.58 (-3.36) 


-0.f4 (-O.fS) 


04 


-0.06 (-0.03) 


-4.99 (-4.69) 


0.30 (0.f5) 


05/06 


0.07 (-0.02) 


-2.96 (-2.57) 


-0.20 (0.06) 


05 


-0.f7 (-O.Of) 


-2.07 (-2.f3) 


0.35 (0.02) 


c" (C/m^) 


el (C/m^) 


gW«rect (C/m2) 


g^rect (C/m^) 


e" (C/m^) 


el (C/m^) 


indirect (C/m^) 


c*— ' (C/m^) 


0.i7 


-S.Of 


-2.85 


-2.58 


-O.fO (0.07) 


4.09 (0.4f) 


3.99 (0.45) 


4.08 (0.48) 



TABLE II: The phonon modes that contribute significantly 
(in the table the percentage contribution is denoted by C) 
in the piezoelectric response. Frequencies are listed for the 
equilibrium geometry and that under -|-f% and -f % strain. 



LAO/PTO 


KSO/PTO 




a; in cm ^ 




uj in cm ^ 


C 


-f% 





+f% 


C 


-f% 





+f% 


f7 


37f.73 


363.04 


352.4f 


fO 


286.24 


284.69 


268.8f 


fO 


3f7.f9 


302.56 


295.93 


29 


2f0.94 


209.94 


2f2.97 


28 


f85.79 


f8f.f8 


f79.5f 


24 


200.6f 


f94.37 


f96.99 


35 


ff6.57 


ff5.f9 


ffO.25 


28 


ff7.44 


f25.83 


ff5.37 



individual superlattices, with the exception of La and 
02/03. 

Due to non-switchable polarization, it is expected that 
these superlattices would exhibit piezoelectric properties 
that show weaker temperature dependence than the con- 
ventional, bulk piezoelectric materials like PbTiOa. To 
explore this, we determine F-point phonons for the su- 
perlattices exhibiting strong piezoelectric responses, i.e 
LAO/PTO and KSO/PTO, and decompose the displace- 
ments of the ions upon strain in terms of these nor- 



mal modes. Table II lists the frequencies of phonon 
modes that dominantly contribute to the piezoelectric 
coefficients and their changes upon +1% and -1% lattice 
strain. The frequencies are found to change by 1-4%, 
except for the softest ones for which the change about 
6-8%. This is in contrast to larger ( 14% ) change in the 
soft mode (u; = 141 cm^^) of tetragonal PbTiOa with 
1% strain, which is known to cause strong temperature 
dependence in its piezo and dielectric responses. With 
relatively smaller Griineisen parameter, these superlat- 
tices, therefore, should have a piezoelectric response that 
is relatively robust as a function of temperature. 

In conclusion, we presented the analysis of piezoelectric 
properties of bi-component, ultra-thin superlattices pro- 
cessing polar discontinuities. We demonstrated tunabil- 
ity of their piezoelectric response via appropriate choice 
of different components, providing a chemical control. 
Moreover their piezoelectric properties are expected to be 
only weakly temperature dependent compared to conven- 
tional piezoelectrics, opening up a new avenue for piezo- 
electrics by design. 
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